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Interesting spin Hamiltonians can be realized with ultracold atoms in a two-component Mott 
insulator (2CMI) [HE]. It was recently demonstrated that the application of a magnetic field 
gradient to the 2CMI enables new techniques of thermometry [3] and adiabatic cooling [4]. Here we 
present a theoretical description which provides quantitative analysis of these two new techniques. 
We show that adiabatic reduction of the field gradient is capable of cooling below the Curie or Neel 
temperature of certain spin ordered phases. 

PACS numbers: 05.30.Jp, 37.10. Jk, 37.10.De, 75.30.Sg, 75.10.Jm, 07.20.Dt 



The possibility of using ultracold lattice-trapped gases 
as general simulators of strongly interacting many-body 
systems has excited increasing interest in recent years [TJ 
[21 [5J. Spin hamiltonians are a natural candidate for quan- 
tum simulation, especially given the relevance of doped 
antifcrromagnetic systems to the important open prob- 
lem of high-T c superconductivity 6 . The 2CMI is the 
starting point for simulation of electronic spin systems 
in lattices [TKTTj. Spin-exchange-stabilized magnetically 
ordered states are expected to exist in the 2CMI [12], 
and observation of these states and transitions between 
them would open up an exciting new field at the inter- 
section of atomic and condensed matter physics. The 
main obstacle which has so far prevented the observation 
of spin-ordered states in the 2CMI is the very low tem- 
perature scale required for spin ordering [13 . Quantum 
Monte Carlo calculations have predicted Curie and Neel 
temperatures on the order of 200 pK for the ferromag- 
netic and antiferromagnetic states of 87 Rb in a 532-nm 
lattice [14 . This is a lower temperature than has ever 
been measured in any system. Clearly, new methods of 
thermometry and refrigeration are required. 

The recently demonstrated technique of spin gradient 
thermometry [3] should allow measurement of tempera- 
tures down to the spin exchange scale in the 2CMI. The 
related method of spin gradient demagnetization cooling 
is capable of cooling to the neighborhood of the critical 
temperature for spin ordering [3]. Together, these new 
techniques open a realistic prospect of preparing spin- 
ordered phases in the 2CMI. In order to compare exper- 
imental results with theory, we have developed a simple 
theoretical model of the 2CMI and used it to calculate 
the expected response of our system to spin gradient ther- 
mometry and spin gradient demagnetization cooling. 

Our treatment of the 2CMI is similar in approach to 
the studies of cooling in the one-component Mott insu- 
lator presented in Refs. |T3] and [15] . in that it is based 
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on calculations of entropy-versus-temperature curves for 
various values of control parameters. Our model neglects 
the effects of tunneling and treats each lattice site sep- 
arately, yet is capable of qualitatively reproducing ob- 
served cooling curves using only one fit parameter (the 
initial temperature) [4]. Our results thus complement, 
and are in qualitative agreement with, the classical mean 
field and Monte Carlo analysis of Natu and Mueller [TB] . 

The inputs to the calculation are the measured trap 
frequencies (u> x , uj y , uj z ) , the total atom number TV, and 
the applied magnetic field gradient V|B| (along with 
various fixed parameters like the scattering lengths and 
magnetic moment of the atoms and the lattice con- 
stant). This allows direct comparison with experiment. 
The particular trap frequencies assumed here are 27r x 
(40, 156, 141) Hz. We assume an atom number of 17000, 
leading to an occupation number of 3 in the center of the 
cloud. These values were chosen because they are typical 
in our experiments. The scattering lengths we assumed 
are Off = 100.4 ao, a^ = 98.98 ao, and = 98.98 ao, 
where ao is the Bohr radius and states f an d i arc 
the \F = l,mp = —1) and \F = 2, rap = —2) hypcrfme 
states of 87 Rb; these values represent the results of the 
most recent theoretical calculations available [17] . 

Detailed technical descriptions of spin gradient ther- 
mometry and spin gradient demagnetization cooling are 
presented in Refs. [3] and [4], respectively. Both tech- 
niques are based on the 2CMI in a magnetic field gradi- 
ent. Since the two components have different magnetic 
moments, the gradient pulls them towards opposite sides 
of the trap, creating two spin domains which remain in 
thermal contact. At zero temperature, there will be a 
zero-width boundary between the two domains, but at 
finite temperature a mixed region composed of spin ex- 
citations will be present. 

Since the total magnetization is always chosen to be 
zero, the average value of the magnetic field is cancelled 
by a Lagrange multiplier and can be subtracted from 
the real field B(x). This allows us to write the field as 
B c ff — V|B| • Xi, where Xi is the vector from the trap 
center to lattice site i projected along the direction of 
the magnetic field gradient. Note that B c g = at the 
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trap center. If tunneling is neglected, then at a magnetic 
field gradient V|B| the energy of a configuration with 
up spins and down spins at lattice site i is 

Ei(n^,n u V|B|) = p ■ V|B| • Xi ■ (n t - n x ) 

cr 

+ • (n T + nx) - A*t"t _ H n ii (!) 

where p is the amplitude of the effective magnetic mo- 
ment of the atoms, Xi — \xi\, a = {til}, U a b is the 
interaction energy between spin a and spin b, Vi — 
(m/2)(cj^xf +UyUi -\-u\zf) is the optical trapping poten- 
tial at site i, yi and Zi are the distances of site i from the 
trap center in the two directions transverse to the gradi- 
ent, and fj, a is the chemical potential of spin a. The chem- 
ical potential is set by the requirement that the number 
of atoms of each spin be equal to half the total experi- 
mentally measured number. 

Equation [l] can be used in the grand canonical 
ensemble to infer the thermal probability of differ- 
ent occupation numbers of the two spins. The 
partition function at lattice site i is Zi(V|B|) = 

E{„ T ,^} ex P(~^K' n i' V l B D/ fc s T )' where k B is 
Boltzmann's constant, T is the temperature, and the 
summation is over all possible combinations of n-j- and 
(each combination is counted only once, due to indis- 
tinguishability of the atoms). The probability of having 
n-f up spins and down spins at lattice site i is then 

Pi (n,, H ,V\B\,T) = e*P(-^K,y,VlB])/ fca T) (2) 

and the resulting entropy at site i is 

$(V|B|,T)= ]T -ftlogp, (3) 

{™T>"4.} 

where the summation is performed in the same way as for 
the partition function. The only additional approxima- 
tion needed is a truncation of the sums over spin configu- 
rations. For our experimental parameters, configurations 
corresponding to a total atom number per site n greater 
than 4 can be neglected, and we have truncated the sums 
accordingly. This truncation is reminiscent of, but more 
general than, the particle-hole approximation |13[ I15j . 
The site entropy Si of Eq. [3] is summed over all lattice 
sites to extract the total entropy as a function of temper- 
ature and field gradient. From this output one can ex- 
tract column-integrated images (Fig[T]), entropy- versus- 
temperature curves (Fig. [2]), and the predicted response 
to thermometry (Fig. [3]) and cooling (Fig. |4| . 

It is instructive to compare the results of this calcu- 
lation to those of the simple approximation which treats 
the spin and particle-hole degrees of freedom separately. 
In this approximation, the partition function for an indi- 
vidual lattice site i is assumed to factorize as Z — Z a Z , 
where Z a = ^ exp(-/9p CT -B(a; i )), is l/k B T, p a is the 




FIG. 1: Comparison of simulated and measured spin images. 
Simulated images are on the left. Magnetic field gradients 
and temperatures for simulated images are: a: 0.7 G/cm, 6 
nK, b: 0.06 G/cm, 2 nK, and c: 0.0024 G/cm, 0.4 nK. The 
gradient and fitted temperature for each measured spin image 
d-f are similar to the values for the simulated image in the 
same row. See Fig. [3] for a comparison of the temperature 
extracted from this fit to the real modeled temperature. Note 
that the total magnetization in all pictures is close to zero. 



magnetic moment of the spin cr, and Zq is the parti- 
tion function of the particle-hole degrees of freedom (for 
which see [T31 HI])- This simple treatment is valid for 
the case of one atom per lattice site. For occupation 
number n > 1, there are corrections which are captured 
by our more complete model. The first correction arises 
from a difference AU between the mean of the intra- 
spin interaction energies U a and the inter-spin interac- 
tion energy f/fj.. AU/U a is about 0.007 for our states. 
The leading correction changes the magnetic field gradi- 
ent at the center of the sample B' to an effective gra- 
dient B' (\ + (n — l)A[//fcsT). This becomes important 
at low temperatures, and destroys the factorizability of 
the partition function mentioned above. The second cor- 
rection is due to indistinguishability of the atoms. This 
arises from the quantum mechanical fact that there are 
three (rather than four) possible spin states for a lattice 
site with two pseudospin-1/2 atoms. The size of both 
corrections is expected to be small, but in order to treat 
them fully we have developed the more general model 
described above. 

Under the assumption that Z = Z a Z§, the mean spin 
(s) as a function of position, field gradient, and temper- 
ature has the simple form 

(3) =tanh(-/3- Ap-B (2:0/2), (4) 

where Ap is the difference between the magnetic mo- 
ments of the two states. Spin gradient thermometry is 
based on the fact that at finite temperatures, the width 
of the boundary layer is proportional to the temperature 
and inversely proportional to the magnetic field gradient. 

The spin profile of Eq. [4] is exact for a 2CMI with 
one particle per site. The model presented here can be 
used to investigate corrections to thermometry at higher 
filling. Figure [3] shows the temperature measured by fit- 
ting to Eq. [4] divided by the actual modeled temperature 
for two values of the gradient. The high-temperature 
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FIG. 2: Total entropy per particle versus temperature, at 
various gradients, for the experimental parameters described 
in the text. The arrow indicates a possible path followed 
during adiabatic spin gradient demagnetization cooling. 
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FIG. 3: Ratio between fitted temperature and actual temper- 
ature at two different gradients, assuming perfect imaging. 
This shows the effect of corrections due to indistinguishabil- 
ity and unequal scattering lengths. Finite imaging resolution 
will limit the range of temperature that can be measured with 
any given gradient, as discussed in Ref. [3|. 



correction is mainly due to indistinguishability of the 
atoms, and is only important for sites containing 2 or 
more atoms. Note that the fitted spin profile is inte- 
grated along both directions transverse to the gradient, 
so it includes contributions from all occupation number 
domains. Though this correction is conceptually impor- 
tant, it changes the measured temperature by less than 
15% under our experimental conditions. The correction 
at low temperatures is partly due to the fact that the 
scattering lengths U^, U~n, and tZ-j-j, are not all equal. 
This is expected to result in a curvature of the mixed 
region between the two spin domains. This curvature 
arises from a buoyancy effect- the species with greater 
intraspin repulsion will preferentially populate the outer 
regions of the trap. This effect causes a fit to Eq. [4] to 
overestimate the temperature, since a curved boundary 
appears wider after integration along the directions per- 
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FIG. 4: Spin gradient demagnetization cooling. Predicted 
temperature versus final gradient, for several values of the 
total entropy. 




FIG. 5: Entropy distribution during spin gradient demagne- 
tization cooling. The images are slices of the cloud through 
the center. Each pixel represents one lattice site. All plots 
are at a total entropy per particle near 0.3 fcg, and are thus 
representative of the changing entropy distribution during 
isentropic demagnetization. Values of the gradient and tem- 
perature are as follows: a: V£> = 0.5G/cm,T = 3nK, b: 
VB = 0.1G/cm,T= I.5nK, c: VB = 0.02G/cm,T = 0.5nK. 
The ring-shaped structures are the mixed-occupation-number 
regions between Mott domains which carry particle-hole en- 
tropy. Note that these regions are narrower after reduction of 
the gradient. This indicates that entropy has been transferred 
from the mixed-occupation-number regions to the mixed-spin 
region, and the temperature has been reduced. 



pendicular to the gradient. Another correction at low 
temperatures arises if the width of the mixed region be- 
comes much less than one lattice constant. In this case 
both the model and the real physical spin system will 
not respond measurably to small changes in the gradi- 
ent, and the measured temperature will overestimate the 
real temperature. These corrections need to be taken 
into account for precision temperature measurements at 
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extreme temperatures and field gradients, but they do 
not alter the conclusions of Refs. [3] or [4]. 

Spin gradient demagnetization cooling is based on the 
fact that the entropy stored in the mixed region between 
the two spin domains increases with decreasing magnetic 
field gradient. If the change of the gradient is adia- 
batic, then the energy and entropy which flow to the 
spin degrees of freedom must come from other degrees 
of freedom, and the sample's temperature can be re- 
duced. Although spin gradient demagnetization cooling 
was inspired by (and is locally similar to) adiabatic de- 
magnetization refrigeration in condensed or gaseous sys- 
tems [TM20] , there are important differences between the 
techniques. Most notably, spin gradient demagnetization 
cooling varies a magnetic field gradient rather than a ho- 
mogeneous field, and relies on spin transport rather than 
spin flips. These differences allow the technique to be ap- 
plied to lattice-trapped ultracold atomic systems. Spin 
gradient demagnetization cooling thus broadens and ex- 
tends existing magnetic refrigeration techniques. 

Entropy versus temperature curves such as those plot- 
ted in Fig. [2] can be used to calculate the response of 
the system to spin gradient demagnetization cooling. If 
the gradient is reduced perfectly adiabatically, the sys- 
tem will move horizontally as indicated by the arrow in 
Fig. [2j and the temperature will decrease. This behavior 
is plotted in Fig. [4] for several values of the total entropy 
(corresponding to different initial temperatures). These 
predictions can be used directly to fit experimental data, 
with the initial temperature being the only free parame- 
ter. Such fitting gives reasonable agreement (see Ref. [1]). 

For sufficiently low initial entropy, the spin degrees of 
freedom will contain all the entropy in the system when 
the gradient is adiabatically reduced by some factor. Fur- 
ther reduction of the gradient below this point is expected 
to linearly decrease the temperature of the system until 
the point where interactions become important. Con- 
versely, if the initial entropy is too high, the spins will 
become fully disordered at some finite value of the gra- 
dient and will no longer be able to absorb entropy. Re- 
duction of the gradient below this point will not change 
the temperature. This behavior, which is essentially a 
finite-size effect, is apparent in the upper curve in Fig. [4] 
If the gradient is sufficiently high, it can pull the two spin 
domains so far apart that the area where they overlap is 
decreased in size. This effect reduces the entropy capac- 
ity of the spin degrees of freedom at high gradients, and 
is the origin of the slight downturn in temperatures at 
the highest gradients in Fig. [4] 

Magnetic field gradients of 1 mG/cm are well within 
the range of the experimentally achievable. Assuming re- 
duction of the gradient from 2 G/cm to 1 mG/cm, our 



analysis predicts that samples with an initial entropy 
lower than about OAks can be cooled below the spin 
ordering temperature. Our model neglects spin correla- 
tions, so the lowest-temperature results plotted in Fig. [4] 
should be taken as evidence that reduction of the gradient 
is capable of cooling below the spin ordering temperature 
rather than a prediction of the dependence of tempera- 
ture on gradient below the Curie or Neel temperature. 

Figure [5] shows several images of the total entropy dis- 
tribution at different final gradients during demagnetiza- 
tion. The pumping of entropy from the kinetic degrees 
of freedom to the spins is clearly visible. 

These theoretical results help elucidate some limita- 
tions on and possible extensions to the technique of spin 
gradient demagnetization cooling. The technique clearly 
requires the use of two states with different magnetic 
moments- this excludes, for example, the two lowest 
states of 7 Li at very high fields. The predicted behav- 
ior is also in principle different for higher filling factors 
than it is for n — 1 (although as discussed above we find 
this effect to be small for the particular case of 87 Rb). 
For example, strong miscibility or immiscibility of the 
two species would change the response of the system to 
demagnetization, but only if the maximum occupation 
number n is greater than 1 (see also Ref. [21] ). The 
dependence of the response to demagnetization on the 
trap frequencies and total atom number can also be in- 
vestigated using the techniques presented here; the most 
important effect of varying these parameters is gener- 
ally to change the maximum occupation number and the 
spectrum of particle-hole excitations. For best cooling 
performance, the initial entropy should be lower than 
the maximum mixing entropy (fcslog2 per site for the 
n = 1 case). We believe that spin gradient demagneti- 
zation cooling could in principle be applied to fermionic 
mixtures as well. In fact, the technique does not even 
require a lattice, and could potentially be applied to the 
thermal fraction of a trapped two-component gas. 

We have presented results of calculations based on a 
theoretical model of the 2CMI and its response to a vary- 
ing magnetic field gradient. Our results provide quanti- 
tative support for spin gradient thermometry and spin 
gradient demagnetization cooling. 
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